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We use the fact that some linear Hamiltonian systems can be considered as 
"finite level" quantum systems, and the description of quantum mechanics in terms 
of probabilities, to associate probability distributions with this particular class of 
linear Hamiltonian systems. 

1 Introduction 

(Sj . Some general structural features of dynamical systems can be exhibited in a study 

of the relation between classical and quantum mechanics. The most widely accepted 
formal relationship between quantum mechanics and classical mechanics is expressed in 
the analogy between commutator brackets of operators on some Hilbert space and Poisson 
brackets among functions on some symplectic phase space, along with the analogy between 
Heisenberg equations of motion and Hamilton equations of motion. This analogy can be 
made more suggestive if operators are replaced by expectation values (quadratic functions 
on the Hilbert space) and the commutator is replaced by Poisson brackets among these 
quadratic functions It is also possible to show that Schrodinger equations define 
Hamiltonian systems on the Hilbert space, where the imaginary part of the scalar product 
defines a symplectic structure. These various aspects have been considered in the past [0 
and more recently in connection with Wigner's problem ||, f|]. They can be summarized 
H ' in the following way. 

On a symplectic vector space, carrier space of a Hamiltonian linear dynamical sys- 
tem, one introduces a complex structure compatible with the symplectic structure, i.e., 
their composition, as maps, defines a positive definite inner product on the vector space. 
Finite-level quantum systems are shown to be linear Hamiltonian systems which, in addi- 
tion, leave invariant a compatible complex structure. Therefore, on a finite-dimensional 
Hilbert space quantum systems are a subclass of linear Hamiltonian systems. This char- 
acterization carries over to infinite dimensions where the "classical system" has an infinite 
numbers of degrees of freedom. By no means, the linear system preserving two compat- 
ible structures (quantum system) can be considered as a quantization of the associated 
Hamiltonian one (classical system). 

Continuing in the spirit of the approach by Moyal ||, Wigner Husimi 0, et al., 
according to which quantum mechanics can be formulated in a natural manner in terms 
of functions on the classical phase space, recently || a different approach has been pro- 
posed dealing with probability distributions rather than quasi-probabilities of previous 
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approaches. It is to be remarked that this approach differs from the one initiated by 
Koopman |§, who had shown how the dynamical transformations of classical mechan- 
ics considered as measure preserving transformations of the phase space, induce unitary 
transformations on the Hilbert space of functions which are square integrable with respect 
to a Liouville measure over the phase space. It is to be stressed that this Hilbert space 
corresponds not to the space of state vectors in quantum mechanics but to the Hilbert 
space of operators on state vectors (with the trace of the product of two operators being 
chosen as the scalar product). 

In this paper, we would like to consider classical canonical linear systems identified 
with the "Hamiltonian part" of quantum systems with a finite number of states. We use 
previous correspondence with probabilities, to associate probabilities also with classical 
linear systems. Then we move to infinite dimensions and show how things appear there. 
This study shows the main difference between classical systems and quantum systems at 
the level of associated probabilities. 

In a future paper, we shall consider how the dynamical evolution and dynamical 
symmetry groups are represented on these probability distributions. 

The state of a classical particle (one degree of freedom) is described by two c-number 
valued functions, position q and momentum p |TD[. In classical statistical mechanics, 
the state of a particle is described by the probability distribution on the particle's phase 
space [JITJ f(q,p) > such that 

J f(q,p)dqdp= 1. (1) 

In quantum mechanics, pure states of a particle are described by a complex wave function 
*&(q) and mixed normalized states are described by density matrices, p(q, q') is the matrix 
element (q | p | q 1 ) of an Hermitian nonnegative density operator p with Tr p = 1. The 
diagonal elements of the density matrix p(q, q) determine the probability distribution in 
position which for pure state p^ —\ | is reduced to p<&(q,q) = \^(q)\ 2 . The wave 

function *&(q) may be considered as a component \l/ g = *&(q) of the state vector in Hilbert 
space. The Hilbert space structure is usually considered as appropriate for description of 
the quantum motion. 

On the other hand, Koopman || has demonstrated that a Hilbert space structure can 



be associated to a classical system as well. This approach has been developed in [[FT 
Also Strocchi || has shown that quantum motion (Schrodinger equation ]12| for the 
wave function) may be reformulated as formally classical motion of coupled harmonic 
oscillators. The linear dynamics of such oscillators and its possible deformations were 
discussed in []T3J, [| |T|. 

The relation of quantum dynamics to a probability distribution evolution (probability 
representation of quantum mechanics) was found in |TJ| in connection with tomography 



transformation [15|, |I^j (Radon transform) of the density matrix used for measuring quan- 



tum states fl7|, [18|. In probability representation, a quantum state is described by a family 
of positive probability distributions of position instead of wave function or density matrix. 
Due to this the probability representation can be considered as classical-like formulation 
of quantum mechanics. The connection of irreducible representations of Heisenberg-Weyl 



group with the tomography transformation was discussed in [19|. 

The goal of our paper is to clarify some mutual relations between these different 
aspects of classical and quantum dynamics. We demonstrate that the classical motion of a 
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system, both with finite and infinite numbers of degrees of freedom, can be associated with 
the probability measure (a family of marginal distributions) characteristic for quantum 
mechanics. 

The paper is organized as follows. In Section 2, we review description of quantum 
dynamics in terms of classical formalism. In Section 3, we consider two- level system in 
detail. Section 4 is devoted to treating two-level system as a spin-l/2-particle system. 
In Section 5, the probability representation (tomographic map) of quantum mechanics 
of spin-1/2 particle is reviewed. In Section 6, the map of states of two and N classical 
oscillators onto probability distributions is considered. In Section 7, we discuss classical 
states of infinite number of oscillators in terms of the marginal distributions and properties 
of the probability distributions under action of canonical transform of the oscillators' phase 
space. 

2 Classical Formalism for Quantum Dynamics 

The Schrodinger evolution equation (h = 1) for the wave function of a one-dimensional 
particle 

i4f(q,t) = m(q,t) , (2) 

where the Hamiltonian operator is 

H = ^ + V(q), P=-<^, Q = Q, (3) 
in the coordinate representation for real and imaginary parts of the wave function 

*(q,t) = %=-j=(x q + ip q )) **(q,t) = ** 9 = -j=(x q -ip q ) (4) 
can be rewritten in the form || 

% ( ) = ( " -H* qq , )(**,)' (5) 
We treat H qq i as matrix elements of the Hamiltonian in the coordinate representation 

(q | H | q') = H qq i (6) 
and perform summation (integration) over repeated indices q' . By introducing notation 

& = (7) 

we rewrite the Schrodinger equation in the form 

^ = hqq'iq' , h qq > = ^ J 9 ^ . (8) 

The evolution equation can be obtained as a classical equation using the Hamiltonian 

n = %H qq ,y q , (9) 
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with summation (integration) over repeated indices and with Poisson brackets 



{%, MV} = {% %} = ; = i5 qq/ . (10) 

The symbol 5 qq > is either Kronecker symbol (if q and q' are discrete indices) or Dirac 
delta- function (for continuous indices q and q'). Thus, Schrodinger equation is equivalent 
to the Hamiltonian "classical equation" 

* fl = {W,¥J. (11) 

One can see that Schrodinger equations are equivalent to classical equations for a system 
of coupled oscillators, with coupling of the oscillators being described by matrix elements 
H qq i. To show this more explicitly, we rewrite the equations for real variables x q ,p q , for 
which Poisson brackets are real ones 

{x q ,x q r} = {p q ,p q >} = 0; {p q ,x q/ } = i5 qq > . (12) 

Then the Hamiltonian takes the form 

1 

-i 

2 

where 



X 9 J 2 \ M - fgg' J ^9?' + H* qq , 

Then the evolution equation for the real vector Q q reads 

Qq — A qq iQ q i , (15) 

with summation (integration) over q' and the matrix A qq i 

A qql = -EB qql ; E = ( _ : . (16) 

We see that the search for the energy levels of the quantum system is equivalent to the 
search for normal modes and their frequencies for a corresponding classical system of 
coupled oscillators. 



3 Two-Level System 

Till now we used infinite system of coupled oscillators and our manipulations were 
done at a formal level. To clarify the approach for constructing probability distribution 
associated to a classical state, let us consider a two-level system. In this case, following [[| 
we have 

* = (£) md 5= (**)- (17) 

The Hermitian matrix H such that 
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is the Hamiltonian and Schrodinger equation looks as 

ii = h = ( % °_ ^ ^ . (19) 

Energy levels of a two- level system are determined by matrix elements a,b = 61 + ib^c of 
matrix (|18|) and we arrive at 



E 1 = £±£ + 1^/(0 + 0)2 + 4(66*-^) 
a + c 1 



E 2 = — -^J (a + c) 2 + 4(bb* - ac) 

By introducing the real components of the wave function 



i.e., 



and 



Q= [1 , (22) 



we arrive at the Hamiltonian description of two formally classical coupled oscillators in 
the form 



1 4 



W = - X! Qa.B a pQp , (23) 

",/3=l 



where 4x4- matrix reads 



1/ H + H* i(H*-H)\ 

2\-i(H-H*) H + H* J' K > 

The Schrodinger evolution equation for the two-level system reads 

Q = AQ, (25) 

where 

A = -EB. 

Thus, the quantum linear dynamics of two-level system is equivalent to the classical 
dynamics of two coupled classical oscillators. This means that one can associate to two 
coupled classical oscillators the two-dimensional Hilbert space of states which formally is 
equivalent to the artificial spin-1/2 states. 

The classical states of two oscillators is described by four real numbers qi,pi, ?2)P2- 
The Hilbert structure is introduced by relating these four numbers to "spinors" 



(26) 



On the other hand, the distribution functions which can determine completely the spinors 
were considered in pO| . Thus, we can associate with two linear classical coupled oscillators 
a positive probability distributions (see below). 



5 



4 Two-Level System as Spin- 1/2 System 



Let us describe now a two-level quantum system and its states in terms of spinorial 
states. To proceed in this direction, we construct the spinor \1/ as a column consisting of 
two-state vectors (of the two- level system) ^ and ^2 as given by formula fllT]). Then the 
density matrix of the two-level system can be constructed for the pure states as 

p<*> = i,k = l, 2. (27) 

For mixed state, the density operator has the form 

P = W 1 \ I +W 2 I tt 2 )<¥ 2 I, (28) 

where W\ and w<i are positive probabilities and u>i+w 2 = 1. The generic Hermitian density 
matrix of the mixed state of the spin-1/2 

\ P12 P22 / 

has the property pn + p 22 = 1- 
For pure state, we have 

pL^I^I 2 ; pS^I^I 2 ; pS? = *i* 2 (30) 

Thus, we can associate to classical two-mode oscillator and to its states (with positions 
gi,g 2 and momenta pi,p 2 ) the density matrix of spin-1/2 particle. 



5 Tomographic map for spinors 



Tomography of spin states was considered in EQ, 21]. The marginal distribution which 



determines the spin state can be constructed using the relation 

Pram' (<P, 9, V) = (DpT>^) mrn , . (31) 

Here p is a density matrix in the initial reference frame and p mm '(<p, 9) is the density matrix 
in the rotated reference frame with Euler angles (<p, 9, ip). The matrix T> is the matrix of 
the spinor irreducible representation of the group 577(2) and its matrix elements depend 
on the Euler angles ((p,9,ip). In view of the structure of the formula for rotated density 
matrix, the diagonal elements of the density matrix depend on two angles (<p, 9) only 
determining the direction of the quantization axis. The one-dimensional stability group 
of p in SU (2) determines the quotient S 2 = SU (2)/U(l). The spin projection on this axis 
has the probability distribution function 

w(m,<p,8) = p mm (Lp,9). (32) 

The function w (m, if, 9) is positive and normalized 

1/2 

w(m,(f,9) = l. 

m=-l/2 
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Since the probability distribution function depends on two extra parameters (<p,9), it is 
possible to express both diagonal and nondiagonal elements of the initial density matrix 
which determines the spin quantum state in terms of the distribution function : 



X 



ir' 2 E e (2j3 + i) 2 

1/2 

E y(-i) m M™i.M)tfS&fa0,VO 



mi =-1/2 



where 



J J JsW J 3 33 \ ^_ _ j.a 
mi —rrii / V m\ —m' 2 m 3 / Sn 2 

-2-K fK r2n 



P%^> ( 33 ) 



o 



du= dip sin 9d9 / dip . (34) 



In formula (|33D, Wigner 3j-symbols and matrix elements of the irreducible representation 
of the SU(2) group Pom are adopted from [IIJ. Thus, the states of spin-1/2 particle are 
mapped onto set of marginal probability distributios w (m, <p, 6) . 



6 Map of Classical States of Two-Mode Harmonic 
Oscillators onto Probability Distribution Set 

In the previous sections, we showed that the states of a classical two-mode oscillator 
labelled by four real c-number functions qi,q2,Pi,P2 are mapped onto state vectors of 
two-dimensional Hilbert space of spinors 

Ql,l2,Pl,P2 < — ► * • 

Spinors, which differ by a common phase factor, i.e., ^ and e lip ^, determine the same 
quantum state. Thus, we arrive at description of pure quantum state of spin-1/2 in terms 
of nonnegative Hermitian 2x2 density matrix p with trace equal to unity and condition 

Tr p 2 = 1 . 

This means that states of classical two-mode oscillators are mapped onto a set of density 
matrices of the pure spin states, 



Qi,Q2,Pi,P2 -»p. 



In Section 5, we have shown that the density matrix of spin state determines the proba- 
bility distribution function 

p — > w (m, (p, 9) 

of the spin projection m onto the axis with direction given by angles ip and 9. Also, 
given the distribution function w (m, ip, 9) one determines the density matrix p due to 
formula : 

w (m, <p, 9) — > p . 
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Thus, we have proved that the classical state of a two-mode oscillator, labelled by four 
real numbers <?i , <?2 , Pi , P2 can be associated with probability distribution, 

Qi,Q2,Pi,P2 -»■ w (m,tp,0) . 

For the spin j, the construction is similar and formulae for spin- 1/2 case are modified as 
follows 



V *2;+i J 

The density matrix of pure state is given by 

i,k = l,2,...,2j + l. 
The density operator of a mixed state has the form 



JJ) 
Pik 



viz -VEr* 



P U) = E I I, 



2j+l 



(35) 



(36) 



(37) 



fc=i 



where 



u k > and 



2j+l 
fc=l 



The rotated density matrix is given by the relation 



(38) 



Here matrices T>^ are the known matrices of irreducible representation of the group 577(2) 
with spin j. For the case of j 7^ 1/2, the formula which connects the density matrix with 
marginal distribution has the natural generalized form of the spin- 1/2 case: 



2i 



33 



(-i)^E E (2j3 + i) 2 

j 3 =0 m 3 = -j3 



X 



X 



E _/(-irMroi,¥>,0) D^{ipM) 



mi=-j 

j j J3 

mi —mi 



J 

rn'i 



J3 



-m' 2 m 3 y 87r 



da; 



,Cj) 



(39) 



"1 < ,L 2 

i.e., in all the cases (where j = 1/2 was used) 1/2 has to be replaced by the value of spin 
equal to j. 

Thus, we obtain the generalization of the map of classical states of an N- dimensional 
oscillator onto a set of distributions. We formulate the result: Given N coupled classical 
oscillators whose state is labeled by 27V real number-coordinates on its phase space 

Q = (Qi,Q2, ■ ■ ■ ,qN,Pi,P2, ■ ■ ■ ,Pn) ■ 
Let us construct complex TV-vector 

( *i \ 
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where 



= —7= (qk + iPk) ; 



k = l,...,N. 



Let us consider now this vector as basis of 2j + 1 = iV-dimensional irreducible represen- 
tation of the group SU(2) with spin j — (N — l)/2. Then points Q in the classical phase 
space are mapped onto the set of probability distribution functions, 



where m = —j, — j + 1, . . . , j — 1, j and angles if, 9 on S 2 . The map is invertible up to 
a common phase factor of the vector \1/ which means we map one-parameter family of 
classical states of 2iV-dimensional coupled oscillators 



with common phase, describing a common rotation in phase space of all the constituent 
oscillators, onto the set of probability distribution functions. 

For iV coupled oscillators, one can consider the limit N — > oo, which implies limit 
j -> oo. 

For spin j, projections m = —j, . . . , j mean that the energy spectrum which is mim- 
icked by frequencies of the normal modes is an equidistant one. But this restriction may 
be easily removed. 



7 Map of Infinite-Dimensional Coupled Oscillator onto 
Distribution Set and Canonical Transformations 



We consider now the system of coupled oscillators which is obtained as a limit iV — > oo 
(j — > oo) of a finite number of coupled oscillators. To do this, we notice that to achieve the 
case of nonequidistant frequency spectrum for "coupled classical oscillators" corresponding 
to nonequidistant energy levels of a quantum system one needs an additional map. This 
map will transform the wave functions of the system with nonequidistant energy levels 
into wave functions of the system with equidistant energy levels. The model of spinors 
is not sensitive to the structure of energy spectrum. The information on the energy-level 
position is an external one and it can be imposed if the physical interpretation of spinor 
components and the corresponding marginal distributions are considered. 

Another aspect to be mentioned is the behavior of the marginal distribution due to 
canonical transformation of position and momentum qk,Pk, i-e., of the real and imaginary 
parts of the wave functions. First, let us consider linear canonical transformations 



where A e Sp (2N, R) is a symplectic real matrix. The Poisson brackets are preserved by 
this transform. If the matrix A has the form 



Q -> w (J) (m,<p,0), 



(Qk + iPk) e 



k + l,2,...,N 



Q' = AQ, 




the symplectic transform can be written as 



(40) 



9 



Here the matrix u and the matrix v are 



u 



Ai + iA 2 A 4 — iA 3 



it 



A4 + ^3 Ai — iX-2 



Formula ( flof) shows the action of the linear canonical transform on wave functions in the 
"classical" formulation of quantum mechanics JJ. 
The density matrix of the pure state 
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is transformed into 



p = u 



(41) 



Formula ( |41~D is an analog of Bogolyubov transform for the density matrix. To show this 
explicitly, we introduce the matrix 

a = ty^ tr . 

Then the transformed density matrix takes the form 

p' = upu + vp*v* + va*u> + uav' . 
For the point canonical transformation defined by 



Xo — A? 







A^ 1 = A 



tr 

4 ) 



one arrives at 



u 



Xi + \\ 



tr 



Ar ltr -A x 



2 ' - 2 ^ 

Since the marginal distribution is related to the density matrix as diagonal element of 
the rotated density matrix, we obtain the formula for canonically transformed probability 
distribution corresponding to the quantum state namely, 



w'(m, <f, 6) 



(43) 



In this way, we have implemented linear canonical transformations on the space of distri- 
bution functions. 



In infinite-dimensional case, we will use the invertible map |22 

,2 



Wy{X,H, V) 



2ttM 



exp 



2u 



iyX 



V(y) dy 



(44) 



of the wave function ^f(y) onto the marginal probability distribution function w^(X, p, v). 
The function w^(X, p, v), with p and v being real parameters, is positive and normalized, 



J Wy(X, p, v) dX = 1 



10 



for normalized wave functions, 



mv)\ 2 dy 



The map can be rewritten for density matrix of pure state p^(y,y') 

pv(v, y') = ^{yW{y') = x y x yl - p y p y < + %{x y p y . - x y >p y ) 

in the form 



W\a(X, p, v) 



2tt v 



P*(y,y') exp 



.y-y 



V 



X-p 



y + y'" 



(45) 



dydy'. (46) 



The Wigner-Moyal function || jj| is related to the generic density matrix by 



W(q,p) 



u u 
P[q+-, q -- 



-ipu 



du 



and 



p(x,x') = —IW 



x + x' 



P e 



ip(x—x') 



dp . 



(47) 



(48) 



On the other hand, in the generic case the marginal distribution is related to the Wigner- 
Moyal function by 



w(X,p,u) = / W(q,p)e 



-ik(X — nq—vp) 



dq dp dv 



and 



W{q, p) = — [w(X, p, v) e W-m-»P) dX dp dv 



(49) 



(50) 



Thus, for pure and mixed states, the probability distribution w (X, p, v) is associated to 
quantum state by means of an invertible integral transform. This property can be used to 
describe quantum states by positive probability distributions instead of wave functions or 
density matrices. In the case of classical pure state labeled by a point in the phase space 
(xg,p g ), the probability distribution w (X,p, v) can be defined, in view of the described 
chain of maps, both for finite and infinite cases. 



8 Conclusion 

We have shown that classical states labeled by points in the classical phase space of 
a classical system can be associated with the marginal probability distribution functions. 
The map constructed uses the map of the points of the classical phase space onto the 
set of wave functions of a corresponding quantum system. Then one uses the invertible 
map from wave functions onto probability distributions introduced in the tomography 
scheme for measuring the quantum states. The established relation clarifies some aspects 
of the connection between quantum and classical pictures of linear dynamical systems. 
This gives the possibility to transport the results of classical considerations into quantum 
pictures and vice versa. 
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